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Abstract 

For a single degree of freedom confined mechanical system with given energy, we know that the motion 
is always periodic and action-angle variables are convenient choice as conjugate phase-space variables. We 
construct action-angle coherent states in view to provide a quantization scheme that yields precisely a given 
observed energy spectrum {En} for such a system. This construction is based on a Bayesian approach: 
each family corresponds to a choice of probability distributions such that the classical energy averaged 
with respect to this probability distribution is precisely £"„ up to a constant shift. The formalism is viewed 
as a natural extension of the Bohr-Sommerfeld rule and an alternative to the canonical quantization. In 
particular, it also yields a satisfactory angle operator as a bounded self-adjoint operator 



*Electronic address: jgazeau @ apc.univ-paris7 .frj 
^Electronic address: |kanamoto .rina @ ocha. ac .j p 



1 



I. INTRODUCTION 



Action-angle variables 

Let us consider a a single degree of freedom confined mechanical system described with phase- 
space conjugate variables {q,p). Suppose it conservative. For a given motion its Hamiltonian 
function is fixed to a certain value E of the energy H(q,p) = E. Solving for the momentum 
variable p leads io p = p{q, E). We suppose that we are in presence of a confinement of the 
system where we have periodic motions only. Then two types of periodic motions are possible: 

(i) Libration: phase trajectory is closed, and then q and p are periodic functions of time with a 
same period. 

(ii) Rotation or circulation: phase trajectory is not closed, then p is periodic function of q. 
For either type, introduce the action variable yj] 

J = jp{q,E)dq=J{E), (1) 

where the loop integral is understood as performed over a complete period of libration (resp. 
rotation). This determines by inversion of J = J{E) the function E = E{J). Now we know 
that the Hamilton- Jacobi equation ^ + H (^q, ^) = obeyed by the action S = J Ldt, with 
L = p q — H, has a solution of the type S = W{q, J) — Et (note that the action S should not 
be confused with the action variable J). The time-independent W = W{J,q) = J pdq is the 
Hamilton characteristic function which generates the contact transformation {q,p) ^ (J,!) at 
constant Hamiltonian, where 7 = ^ is the angle variable, conjugate to J. It follows from the 
definition of J that the period r (resp. frequency of motion at fixed energy is r = dJ/dE = 
t{E) (resp. i'{E) = 1/t(E) = dE/dJ), and so the time evolution of the angle variable is linear, 
with period r: 

7 = -4^ + 70 = J^iE) t + 70 • (2) 
Note that this equation allows to consider the time as proper to the system at given energy. 

Action-angle coherent states for measured energies 

Let us suppose that a series of energy measurements on a mechanical system with one-degree 
of freedom yields an energy spectrum Eq, Ei, . . . En, hi this paper families of corresponding 



action-angle coherent states are constructed in view to provide a quantization scheme consistent 
with that discrete sequence of experimental energies. The construction is based on a Bayesian 
approach: each family corresponds to a choice of probability distributions, n i— Pn{J)/.N'{J) 
(prior discrete), J i— )■ Pn{J) (posterior continuous) such that the mean value of the classical energy 
with respect to the probability p„( J) is precisely En, up to the addition of a constant independent of 
n. The formalism \^ can be viewed as a natural extension of the empirical Bohr-Sommerfeld rule, 
J{E) = nh, where h is the Planck constant. We know that this quantization is exact for the motion 
on the circle (quantization of the angular momentum), and valid in the semiclassical regime. In 
the deep quantum regime our approach can be viewed as a viable alternative to the canonical 
quantization, particularly when the latter is impracticable. For instance, it yields a satisfactory 
angle operator. 

In Section HI] we give an overview of a general construction of coherent states labelled by 
elements of a measure set, their use as a quantizer frame, and their suitability in regard with their 
phase space content. We specify this formalism in Section HlI] to the action-angle phase space of a 
confined system by following a Bayesian construction mentioned above. In Section |IV] we revisit 
the question of suitability of these action-angle coherent states. Sections |V] and |VI] are devoted to 
the question of choice of probability distributions appropriate for our quantization goals. Section 
I VIII is an illustration of our approach by examining a few characteristics, like time evolution, of 
two types of coherent states for the free rotator. We end this note in Section IVIIII with some 
comments about questions raised by our approach and possible generalizations. 

II. COHERENT STATES FAMILY AS HILBERTIAN FRAME FOR PHASE SPACE 
A. The general construction iQ] 

Let X be the phase space of a mechanical system, equipped with its symplectic measure /i. 
Actually (X, /i) can be any measure space. Let X be some countable set and O be an orthonormal 
system O = {(pn,n' G X} made of elements 0„(x) in the Hilbert space L^(X, /i), with the 
positiveness and finiteness constraints 



Let Ti be some separable Hilbert space with orthonormal basis {|e„) , n G X} in one-to-one corre- 
spondence with elements of C: |e„) t-)- 0^. Then the following family of states in the companion 




(3) 



n 
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%, labelled by elements of X, 



\^) = ^ffr^y^'^n{x)\en) , (4) 



where </)„(x) refers to the complex conjugate, and obey normalization, {x\x) = 1, and resolution 
of the unity n(dx) J\f{x) \x){x\ = 1^. Due to these two properties, vectors \x) are named 
coherent states (CS) (in a wide sense). The resolution of the unity is precisely the departure point 
for the corresponding CS quantization of functions (or distributions when the latter are properly 
defined) on X which transforms them into linear operators in T-L: 



f{x)t-^Af= / jj,{dx)Af{x) f{x)\x){x\ . 
Jx 



(5) 



f{x) may be considered as CS quantizable if A j is densely defined in H or if the so-called lower 
symbol of Af, f{x) =^ {x\Af\x) is a smooth function on X viewed as a topological phase space. 
Hence, the family x G X of coherent states offers a certain point of view or frame (in 

the right Hilbertian geometrical sense) to analyze, in a non-commutative way, the classical set of 
points X. Changing the family corresponds to a change of frame and possibly to an equivalent or 
non equivalent quantization in regard with specific physical quantities, like energy, action, angle, 
position..., and their mutual Poisson brackets. 



B. What we understand by suitable coherent states in Quantum Mechanics 

Besides the fundamental quantizer role played by the CS further criteria are usually re- 
quested to qualify the latter as suitable from the classical-quantum relation point of view. We list 
here a few of them. For this purpose, we suppose that the CS \x) depend on a parameter, say h, 
the Planck constant, such that the classical limit corresponds to /i — )■ 0. 

(i) Relative error criterion. For a semi -bounded quantizable function f(x) we define the 
relative error: 

f{x)''=' {x\Af\x), (6) 

where C is a constant which is chosen such that |/(a;) + C| 7^ for all x. Hence: a CS 
family is said suitable as a family of quasi-classical states with respect to the function f if 
supj, rerry(a;) as h ^ 0. 



rerr j (x) 



fix) - fix) 



fix) + C 
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(ii) Time evolution criterion in semi-classical regime. Let H{x) be the Hamiltonian of a 
mechanical system with phase space X and let Ah be its CS-quantized counterpart. We 
define the time evolution of the probability density of a CS \xo) in phase space with respect 
to the measure ^{dx) as the function 

p,,{t,x)^^ Af{x)\{x\e-'^'''/'\xo)\\ (7) 

Its probabilistic nature is directly derived from the resolution of the unity. Hence: a CS 
family is said suitable as a family of quasi-classical states with respect to the time evolution 
ruled by the quantum Hamiltonian Ah if the support of the function p^^ (t, x) tends to locate 
on the classical phase space trajectory x(t) with initial condition x{0) = Xq as h ^ 0. 

(iii) Time evolution stability. A CS family is said temporally stable if under the action of the 
evolution operator e.-'^^at/h^ ^ coherent state is transformed into another coherent state in 
the same family, up to possibly a phase factor: 

e-'^^'/''\x) = e'^^'^\x{t)) . (8) 

For more details, see [|4|] and references therein. 



ni. A BAYESIAN PROBABILISTIC [5] CONSTRUCTION OF ACTION-ANGLE COHERENT 
STATES AND RELATED QUANTIZATIONS 

Conditional posterior probability distribution 

Suppose that measurement on a confined one-dimensional system yields the sequence of values 
for the energy observable (up to a constant shift): 

^0 < ^1 < • • ■ < K < ■ ■ ■ • (9) 
Supposing a (prior) uniform distribution on the range of the action variable J, we define a cor- 
responding sequence of probability distributions J ^ Pn{J), i-C- /iRorR+ dJ Pn{J) = 1, with 
J def J ij^^ obeying the two conditions: 

< M{J) = Yl P-^-^^ < °° ' + est = / dJE{J)pniJ) , (10) 

neZorN ^MorR+ 

where R and Z (resp. R"*" and N) stand for the rotation (resp. libration) type of motion. The finite- 
ness condition allows to consider the map n h-> p„( J) /N{J) as a probabilistic model referring to 
the discrete data, which might be viewed in the present context as a prior distribution also. 



Action-angle coherent states 



Let Hhe a complex separable Hilbert space with orthonormal basis {|e„) , n G ZorN}. Ti 
is the space of quantum states. Let r > be a rescaled period of the angle variable and X = 
{( J, 7) , J G M orM+, < 7 < r} be the action-angle phase space for a rotation (resp. libration) 
motion with given energies the discrete sequence Eq < Ei < ■ • ■ < En < ■ ■ ■ . Let iPniJ))nezor-N 
be the sequence of probability distributions associated with these energies. We suppose J) = 
Pn{—J) in the rotation case. One then constructs the family of states in "H for the rotation or 
libration motion as the following continuous map from X into 1-L: 

X 3 (J, 7) ^ I J,7) = Yl VP^e-""" 10 en, (11) 

where the choice of the real sequence n 1— >■ a„ is left to us in order to comply with some if not all 
criteria previously listed. 

Fundamental properties of action-angle coherent states 

In both cases the coherent states | J, 7) 

(i) are unit vector : (J, 7I J, 7) = 1 

(ii) resolve the unity operator in "H with respect a measure "in the Bohr sense" usidJ d'j) [|4[] 
on the phase space X : 



T 

/•-TOO ^ 1 /"^ 

fiB{dJd^)Ar{J)\J,-f){J,^\''^' / dJATiJ) lim - / d7|J,7)(J,7| = lw, 

T-i-oo 1 It 



X 



(12) 

Here, we impose T = 2Mt with M G N and letting M — )■ 00. Hence, if the sequence 
(a„) assumes all its values in 27rZ/r, then the angular integral reduces to ^ d^y and the 
measure fiB becomes the ordinary one on the cylinder. 

(iii) allow a "coherent state quantization" of classical observables /(J, 7), 

/(J, 7)^ [ ^BidJd^)XiJ)fiJ,^)\J,^){J,-f\ = Af, (13) 
Jx 

which is compatible with the energy constraint (flOl) on the posterior distribution J ^ pn{J). 
Indeed, the CS quantized version of the classical Hamiltonian H = E[J) is diagonal in 



the basis {|e„) , n = 0, 1, ... } since it is trivially verified that in both cases the quantum 
Hamiltonian is exactly what we expect: 

^E{j) = J^^En + cst)|e„)(e„| . (14) 

n 

Particular quantizations 

Actually, the quantization of any function factiJ) of the action variable alone yields the diago- 
nal operator: 

/act(^) %ct = X^(/act)n,|e„,)(e„| , (/act)n, = / dJ f^ct{J) Pn{J) ■ (15) 

„ ^EorIR+ 

On the other hand, the quantization of any r periodic function /ang(7) of the angle variable alone 
yields the operator: 

/ang(7) ^ ^/ang = ^l^ujnn'len) (e„' | , (16) 



n,n' 



where the matrix elements are formally given by: 



T 

1 



[AuJnn'= / dJ^pniJ)pn'{J) lira - / rf7e-^(""-""')Vang(7) 

if an-an'^^Z, ^^^^ 

^nn' Cfc(/ang; t) if - «„' = ^A; G , 

where Ck{f]T) = ^ d'y f{^) e~^'^'^^'^^'^ is the kih. Fourier coefficient of /(7), and Wnn' = 
/RorR+ ^/Pn{J) Pn'{J) mcasurcs Correlation between the two distributions J H- Pn{J), J ^ 
Pn'{J)- Note that the diagonal values are all equal to the average of /ang(7) over one period. Also 
the infinite matrix can be sparse, even just diagonal, depending on the choice of the p„( J) and 
a„'s. Hence, the quantization could at end transform classical observables into a commutative 
algebra of operators. 

An important point is that this CS quantization procedure provides, for a given choice of the 
sequence a self-adjoint angle operator corresponding to the angle function ^(7) defined 
on the real line as the r-periodic extension of ^(7) = 7 on the interval [0, r). Then Ck{A] r) = 
i/{2'Kk) for A; 7^ and co(^; r) = r/2. 
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IV. THE QUEST FOR "SUITABLE" COHERENT STATES 



From relative error 

Since the relative error function involves lower symbols, let us just consider those for the two 
particular cases of classical functions, /act(</) and /ang(7)- 

hct{J) = {J,-f\Af^JJ,-f) = X^(/act)n^^|j| = ((/act)n)j, 

n,n' 

where we observe, in the first case, the appearance of a double averaging using the two Bayesian 
facets, and, in the second case, the presence of a deformation of the Fourier series of /ang (7) 
involving, on one hand, the selection rules k{n, n') =^ ^ (a„ — G Z, and weights of proba- 
bilistic origin on the other hand. The relative error ([6]) expresses relative deviations of the original 
classical observables to the above types of averaging involving the two facts of the underlying 
Bayesian duality. 

From localization probability distributions 

The action-angle phase space representation of a particular coherent state | Jq, 70) , as a function 
of (J, 7), is the "normalized" overlap 

*|Jo,7o)(-/,7) = vW) (^,7l^o,7o) = A^ ^Vp„.(J)p„(Jo)e-"(^-^°) , (18) 

Hence, the map X 3 (7,7) ^ p\X"ZM,l) = I^|Jo,7o>(^, 7)1' = ^{J) l(^,7l^o,7o)|' represents 
a localization probability distribution, namely a generalized version of the Husimi distribution, on 
the phase space provided with the pseudo-measure ^b- Indeed, the resolution of the identity gives 
immediately 

[ /i^(dJci7)pfS)(J,7) = l- (19) 

If we choose instead a specific realization of the Hilbert space "H, like that one generated by 
eigenfunctions of the quantum Hamiltonian in "g" or "configuration" representation, |e„) h- >^ 
ipnil), the corresponding representation of the state | J07 7o) reads as 

i^UoaoM = Jrm E VMJo)e-'''"^° Mq) , (20) 



with corresponding probability density of localization on the range of the g-variable given by 
From time evolution 

Since the CS quantized version Ah of the classical Hamiltonian H is diagonal in the basis 
{ |e„) , n E Z (resp. N) }, the time evolution of the CS in both representations is given respectively 
by 

= Yl VPn(JoK(J)e^^""(^-^°^"^"*/'^^ , (21) 

"*/V|Jo,7o)(?) = E V^;:W e-(-^o+^"*/'^) Mq) ■ (22) 



Snapshots of the time evolution of the corresponding probability densities 
|a/A/'(J) ( J, 7|e~*^-f^*/'^| Jo, 7o)P and \e~''^"^^'^tp\jQ^^„){q)\'^ are necessary in order to discrim- 
inate suitable functions Pn{J) and sequences («„) from the classical limit viewpoint. Note that 
temporal evolution stability is granted with the choice a„ = a_„ = En- 

V. THE QUEST FOR EXPLICIT PROBABILITIES n p„ ( J) : THE TWO SIMPLE CASES 

The quantities that are left undetermined in our construction of the CS are the discretely indexed 
probability distribution J 1— )• p„( J) and the sequence n i-t- a„. Two simplest situations are helpful 
in giving some hints: the free rotator (mass m on circle of radius /) and the harmonic oscillator 
(frequency cu), for which the energies are respectively En oc n^+const. and E^ oc n+const. 

In the first case, a familiar solution [|61] is family of Gaussians centered at each integer, with 



dimensionless width parameter ci or e = 1/ {2a ): 

= (2^) - (^) . " ^ Z . (23) 

with J = J/h. This gives the eigenvalues J„ = hn of Aj and 

^ h^n^ ^ a^h^ _ h^n^ ^ ^^^^ 
" Sir'^mP Sn'^mP Sir'^mP IGevr^m/^ ' 

the constant shift being the average value of the classical energy with respect to the distribution 

Po{J). By introducing the Compton length Ac = h/mc of the particle Eq. (|24l) reads E^ = 



{\c/lYmc^{n^ + l/(2e))/2. As an illustration of the suitability of this Gaussian choice, we show 
in Figure [H the respective behaviors of the spectrum and lower symbol of the angle operator 
for different values of the parameter e. 
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FIG. 1: Spectrum (left) and lower symbol (right) of the angle operator obtained by CS quantization of the an- 
gle of rotation on the circle, for different values of parameter e = l/(2cj^), whenp„(J) = (f)^''^ ^-'^i^-'^)'^ 
and a„ = n. For the spectrum, e = 10"'', 0.3, 1, 3, 5, 50. For the lower symbol: e = 10"'', 0.1, 1, 3, 10 
(the lower symbol has little dependence on J). One notices the tendency to the sawtoothed behavior of the 
classical angle function as e — )■ oo, i.e. as cr — > 0. 



In the second case, another familiar solution is the discretely indexed gamma distribution: 

- 7" 

Pn(J)=e-^ — , riGN. (25) 
n\ 

This gives J„ = + 1) and £"„ = ^^^{n + 1) = hwin^X), the constant shift being the average 
value of the classical energy with respect to the distribution po{J). 

VI. THE QUEST FOR EXPLICIT PROBABILITIES : THE GENERAL CASE 

Given the classical relation E = E{J) between action variable and energy, and the obser- 
vational or computed sequence {En), the central question is to find the sequence of probability 
distributions J ^ Pn{J) (at least with a satisfying approximation), which obey the two funda- 
mental conditions (flOl) . For the rotation case, the departure point could be a normal-like law, 
possibly modified along a perturbation scheme with expansion parameter the strength of the po- 
tential energy U. Similarly, in the libration case, the departure could be a gamma-like law, possibly 
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modified along the same lines. A nice and manageable model is the simple pendulum, whose the 
quantum version is well known from the solution of the Mathieu equation Isj]. 



Approximations for the simple pendulum: the rotation case 

Starting from the computed Mathieu eigenvalues En, an empirical approach consists in starting 
from the sequence of computed action variables J^^ =^ J{En) from Eq. ([T]) and to impose, in the 
rotation case, the sequence of normal laws centered at these J^. 

Pni-J)=(7r^] e'^^'-'-\ neZ, (26) 
by "adjusting" (7„ in order to suitably approximate the EnS with the computed quantities 

/ + 00 
E{J)pn{J)dJ. (27) 
-oo 

Note that with this choice, the eigenvalues J„ of Aj axe precisely the J^''s. 



Approximations for the simple pendulum: the libration case 

Handling the libration case is more delicate. Another empirical approach consists in introduc- 
ing the gamma-like distribution: 

1 jn °° jn 

T / T\ L J , dcf I dcf r" / T\ dcf \ ^ 'J 

J ^ Pn{J) = jrjjz , nGN, yj = yiy2 ■ ■ ■ Vn, yoi = l, Sy{J) = 2_^ — 

with an auxiliary sequence {0 = yo < yi < ■ ■ ■ < yn < ■ ■ ■ } such that the corresponding moment 
problem has a solution with a positive measure Wy{J) dJ, 



dJ Wy{J) Pn{J) = 1 , (28) 

and that the quantization conditions involving the computed Mathieu eigenvalues, 

r+oo 

En + cst= dJ Wy{J) E{J) pn{J) , n = 0,l,2,---, 
Jo 

are fulfilled (at least approximately). We note that, by construction, we have J2^=oPniJ) = 1> 



that the action variable J is the mean value of n (-)■ ?/„ with respect to the Poisson-like distribution 

n ^ pn{ 

{J) J = Vn+l- 



n H- Pn{J)' (yn) '= Yl^=oynPn{J) = J, ^ud that the quantum action Aj has eigenvalues 
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VII, FREE ROTOR CS : TWO INTERESTING CHOICES 



With the Gaussian choice (l23l) and for a general choice of a sequence of frequencies an, the 
coherent states for the free rotator read 

|J,7) = (l)'/'^e-t(^'-")%— "^|e„) , (29) 

The normalization function M{J) is given in two forms: 

and satisfies \im^^oJ\f(J) = 1. In the following we consider two choices of a and investigate 
what kind of properties of CS is satisfied for each case. 

Case On = 2-Kn/T 

This choice renders exact the quantization of the classical canonical commutation rule 

{ J, e«27r7/r| ^ j^^i2nj/T _ ludccd, wc havc [Aj, A^i27r-y/r] = hA^i2nj/r. Concerning the phase space 
distribution ^ for time evolution, we obtain the following upper bound: 

p.o.o(^,7;t) < .;^- /XJ7 ^^"'^^'^""''''^^' ^''^ 



where e = 1/(2(t^), fi = {J + Jo)/2, t = ht/{2mP), and 7 = 27r(7 — 7o)/r. From this is derived 

1 

1/1+40-442 



the estimate on the semi-classical behavior at large Jq = M E Z: pjo,yoiJ^ 1'-> ^) ^ / Sjj^ 



for 7/27r - Mt G Z , and vanishes if Jq ^ Z or 7/27r - Mt ^ Z. 
Case a„ = lirn? /t 

This choice is appropriate to temporal evolution stability 

e-'^»*'%J, 7) = e"^*/2 J, 7 - — 1\ . (32) 

27r / 

This time, the phase space distribution ([7]) is bounded as follows. 

A = ^(7-7o-/iV(4vrm/2)) 



Pio.o(^,7;t) < , ' ^T77T^E^"^^^^^"''^'"^ (33) 
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where 7(t) = ^(7 — 70 — ^t/ {AirmP)) = 7 — 1 Hence the estimate on semi-classical behavior 

of Pjo,7o('^. 7; i) at large Jq = M e Z: pjo,^o( J, 7; t) < ^^^/^4^2(t) '^JJo 7(^) e 27rZ/Af , and 
vanishes if Jq ^ Z or 7(t) ^ 27rZ/M. 

Vin. CONCLUDING POINTS 

Periodic and more generally Integrable systems provide a variety of such families of action- 
angle coherent states. A future issue is about the best choice of probability distributions Ui 
Pmi-Ji)- A more fundamental concerns the physical (in terms of physical measurement) equiva- 
lence between different families of coherent states from a quantization point of view. Finally, the 
extension of the probabilistic approach presented in this work to unconfined systems and subse- 
quent continuous spectra is possible ||7|]. 
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